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High resolution spectroscopy of the BCD galaxy Haro 15: 
I. Internal kinematics 
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ABSTRACT 

Using echelle spectroscopy, obtained at Las Campanas Observatory, we present a de- 
tailed study of the internal kinematics of the nebular material in multiple knots of 
the blue compact dwarf galaxy Haro 15. A detailed analysis of the complex emission 
line profiles show the presence of an underlying broad component in almost all knots, 
and the brightest star-forming region shows unmistakable signs for the presence of 
two distinct narrow kinematical components. We also study the information that our 
analysis provides regarding the motion of the individual knots in the Haro 15 galaxy 
potential, confirming that they follow galactic rotation. Finally, we examine the re- 
lation between their velocity dispersion and luminosity, finding that almost all knots 
follow the relation for virialised systems. This holds for the strong narrow components 
identified in complex fits and for single profile fits, although the latter show a flatter 
slope. In agreement with previous findings, in this paper we show that the existence of 
multiple kinematical components among massive starbursts cannot be overlooked, as 
it has a noticeable effect on any subsequent analysis that relies on basic parameters. 
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1 INTRODUCTION 

Already in the original list of blue galaxies with emission 
lines by Haro (1956) who described this object as "minute 
comet ary nebula", Haro 15 has later been included in sev- 
eral compilations of Blue Compact Galaxies (BCG). Several 
studied at all frequencies have been carried out over this ob- 
ject. In particular, the optical spectroscopy was analysed by 
Hunter & Gallagher (1985); Mazzarella et al. (1991); Kong 
et al. (2002); Shi et al. (2005), between others. Owing to the 
detection of the He ii A4686 emission line by Kovo & Con- 
tini (1999), Schaerer et al. (1999) has classified this galaxy 
as a Wolf-Rayet (WR) galaxy. Lopez-Sanchez & Esteban 
(2010) also found a blue WR bump supporting the WR na- 
ture of Haro 15. The Ha image shown in Cairos et al. (2001) 
shows a knotty morphology with the starburst region re- 
solved in a large number of clumps which appear scattered 
over the entire galaxy. The integrated Ha luminosity yields 
a star formation rate (SFR) of 3.3 M© yr~^ (Lopez-Sanchez 
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& Esteban 2008). Lopez-Sanchez (2010) analysed the SFR 
in Haro 15 using mult i- wavelength data. The derived values 
using Ha, FIR (far-infrared) and radio data are very similar, 
but the value derived using the FUV (far- ultraviolet) lumi- 
nosity is twice the others, confirming that a young stellar 
population is dominating the light of the galaxy and suggest- 
ing that the starburst phenomenon in Haro 15 started some 
time ago (at least, 100 Myr ago). From the morphological 
point of view, Haro 15 has been classified as an (R)SBO pe- 
culiar galaxy by de Vaucouleurs et al. (1991) and in fact, in 
the deep images taken with the 2.2m CAHA telescope shown 
in Lopez-Sanchez & Esteban (2008), the spiral morphology 
of the galaxy can easily be appreciated. Two high surface 
brightness concentrations can be distinguished, named A 
and C by the authors. Both show blue colours and high 
FUV emission, indicative of recent and on- going star forma- 
tion activity, further supported by the presence of the WR 
features. 

Long slit optical spectroscopy of Haro 15 has been pre- 
sented in Lopez-Sanchez & Esteban (2009) covering four 
main regions in the galaxy: the center (named C by the 
authors) and three bright regions located ESE (named A), 
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WNW (named D) and NE (named B). The spectrum of 
the central region shows strong nebular emission lines and 
prominent stellar absorption wings in the Hi Balmer lines, 
evidence of an underlying stellar population which the au- 
thors estimate to be around 500 Myr old. Only emission 
lines are seen in regions A and B and the spectrum of re- 
gion D results too noisy to be analysed. These authors de- 
rived an oxygen abundance of 12+log(0/H)= 8.37=b0.10 and 
8.10ib0.06 for regions A and B, respectively. Based on these 
data the authors conclude that all the observed knots can 
be classified as typical Hii regions. 

Haro 15 has an absolute magnitude M^b — -20.69, a sur- 
face brightness fiB— 18.56 mag arcsec"^ and a colour B-V= 
0.33 (Cairos et al. 2001). At a distance of 86.6 Mpc (de Vau- 
couleurs et al. 1991), Haro 15 meets the criteria for a Lu- 
minous Compact Blue Galaxy (LCBG) (Hoyos et al. 2004) 
despite the fact of its brightness distribution showing an ex- 
ponential profile with a scale length of 1.37 kpc. The nature 
of this type of galaxies is not yet clear since they probably 
constitute a mixed population of star burst galaxies. Some 
authors suggest that LCBGs might represent the final out- 
come of a merger between a dwarf elliptical and a gas rich 
dwarf galaxy or Hi cloud (Ostlin 1998; Gumming et al. 2008). 
The interactions taking place during the merging process 
would act as the starburst trigger (see e.g. Lopez-Sanchez 
et al. 2004, 2006). The images of Haro 15 showing two sep- 
arated nuclei surrounded by a more regular, roughly ellipti- 
cal envelope with twisted isophotes, together with different 
features and faint extensions as reported by Gairos et al. 
(2001) give support to this picture. Other galaxy properties 
like different kinematics and chemical abundances between 
regions A and B may also indicate that Haro 15 is probably 
experiencing a minor merger (Lopez-Sanchez 2010). On the 
other hand, a reliable determination of the mass of this kind 
of objects is needed in order to decide on the evolutionary 
path they may follow: becoming dwarf spheroidal galaxies 
(Hoyos et al. 2007) or the bulges of small spirals (Hammer 
et al. 2001). However, interactions and mergers, as well as 
feedback processes result of intense star formation, might re- 
flect in a peculiar gas kinematics that prevent the derivation 
of their dynamical masses from the gas velocity widths. 

Echelle spectroscopy provides a means to look for kine- 
matically different components in the emission lines of the 
ionised gas since it reaches the spectral resolution needed 
to resolve the presence of structures within the emission- 
line profile, usually masked by its large supersonic width. 
In Firpo et al. (2010) we observed a residual present in the 
wings of several lines when fitting a single Gaussian pro- 
files to the emission lines observed in the Giant Hii Regions 
within the galaxies NGG 7479 and NGG 6070. Basing on the 
variety studies that have been proposed in the literature to 
interpret the existence of the broad supersonic component 
measured in the emission line profile of Giant Hii Regions, 
and whenever possible, we have evaluated the possible pres- 
ence of a broad component (Mufioz-Tuhon et al. 1996; Mel- 
nick et al. 1999; Hagele et al. 2007, 2009, 2010, among oth- 
ers) or two symmetric low-intensity components in the fit to 
the observed emission line profile widths (Ghu & Kennicutt 
1994; Relafio k Beckman 2005; Rozas et al. 2006). In the 
present work, we have also found that all Haro 15 regions 
show evidence of wing broadening evident mainly in the Ha 
line and confirmed in other emission lines. Another common 



Figure 1. Finding chart for the Giant Hll Region candidates 
observed in Haro 15, identified by circles and labelled following 
Cairos et al. (2001). Wide Field Planetary Camera 2 Ho: image 
was obtained from the Multimission Archive at the Space Tele- 
scope Science Institute (MAST) 




feature of Giant Hii Regions is the supersonic width of their 
emission line profiles, although the origin of these velocities 
is not yet clear (Terlevich & Melnick 1981; Roy et al. 1986; 
Hippelein 1986). A relation between Ha luminosity and ve- 
locity dispersion of the form L(Hq;) oc cf^ is expected under 
the assumption of a gravitational origin of the ionised gas 
dynamics. Therefore, a study of the implications that the 
existence of several kinematically distinct components has 
on the L(Ha) vs. a relation is of great value. 

In this paper we present echelle spectroscopy obtained 
with the 100-inch du Pont Telescope at Las Gampanas Ob- 
servatory (LGO) of five different positions across Haro 15 
which provide a velocity resolution of about 12 kms~^ allow- 
ing the identification of different kinematical components of 
the gas and the measurement of their corresponding velocity 
dispersions. This study is part of a project to obtain high 
spectral resolution echelle data to determine the nature of 
Giant Hii Regions visible from the southern hemisphere and 
analyse the physical conditions of the ionised gas of these 
regions and Blue Compact Dwarf galaxies. In Section 2 we 
present the observations and the data reduction. Section 3 
presents the analysis of the emission-line profiles and dis- 
cusses the results. Finally, the summary and conclusions of 
this work are given in Section 4. 
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Table 1. Journal of observations. The identifications in the first 
column correspond to the nomenclature used in this paper, fol- 
lowing and extending that in Cairos et al. (2001); columns 2 and 

3 list the observation date and the air mass respectively; columns 

4 and 5 show the distance to the galactic center in arcsecs and 
kpc, respectivey; finally, column 5 lists the identification of the 
regions given in Lopez-Sanchez Esteban (2008) referred to as 
LS notation. The exposure time for each region was 1800 seconds. 



Knots 


Date 


sec z 


Dist.(") 


Dist.(kpc) 


LS 


Haro 15 A 


2006 Jul 19 


1.2 








C 


Haro 15 B 


2006 Jul 19 


1.1 


11 


4.62 


A 


Haro 15 C 


2006 Jul 19 


1.1 


20 


8.4 


B 


Haro 15 E 


2006 Jul 20 


1.2 


16 


6.72 




Haro 15 F 


2006 Jul 20 


1.1 


17 


7.14 





2 OBSERVATIONS AND DATA REDUCTION 

High resolution spectroscopy of Haro 15 was obtained 
using an echelle spectrograph attached to the 100-inch 
du Pont Telescope, Las Campanas Observatory (LCO), 
in July 19 and 20 of 2006. The spectral range covered by 
the observations was from 3400 to 10000 A. Observing 
conditions were good, with an average of 1 arcsec seeing 
and photometric nights. A 2x2 binning was applied to the 
CCD in order to minimise the readout contribution to the 
final spectrum noise. With 1 arcsec effective slit width and 
4 arcsec slit length, the spectral resolution achieved in our 
du Pont Echelle data was R^25000: AA=0.25A at A6OOOA, 
as measured from the FWHM of the ThAr comparison lines 
taken for wavelength calibration purposes. This translates 
into a velocity resolution of ^^12 kms~^. 
Considering the distance to Haro 15 and that the spatial 
resolution is limited by seeing, the smallest structure that 
can be resolved is 0.43 kpc (^ l"'-^ 1.93 pixel) 

Five different regions in Haro 15 were observed as shown 
in Figure 1 . The spectra were obtaind as singles exposures of 
1800 seconds. CALSPEC spectrophotometric standard star 
Feige 110 (Bohlin et al. 2001) was also observed for flux cal- 
ibration purposes with an exposure time of 1200 seconds. In 
addition, Th-Ar comparison spectra, milky flats (sky flats 
obtained with a diffuser, during the afternoon) and bias 
frames were taken every night. A journal of observations 
is shown in Table 1. 

The data analysis was carried out using the IRAF^ 
software. After bias subtraction and flat field corrections 
by means of milky flats the bidimensional images were cor- 
rected for cosmic rays using the task COSMICRAYS which de- 
tects and removes cosmic rays using a flux ratio algorithm. 
The corrected data were reduced by IRAF routines following 
procedures similar to those described in Firpo et al. (2005). 

At the end of the process, we compared the red end of 
the wavelength calibrated spectra with the night-sky spec- 
trum by Osterbrock et al. (1996). This turned out to be a 
very reliable confirmation of the goodness of the wavelength 
solution, and we checked that differences between our wave- 
lengths and the sky line wavelengths were below 0.04 A. 



^ Image Reduction and Analysis Facility, distributed by NOAO, 
operated by AURA, Inc., under agreement with NSF. 



Flux calibration was achieved by observing the CAL- 
SPEC spectrophotometric standard star, Feige 110 (Bohlin 
et al. 2001) whose flux was tabulated every 2A. Despite its 
relatively low brightness (V=11.83), Feige 110 is ideal for 
calibrating high resolution echelle spectra. The amount of 
defined intervals within an echelle order ranged from four 
to twelve, depending on the quality of the spectrum. The 
flux calibrations was performed as described in Firpo et al. 
(2005). 



3 RESULTS AND DISCUSSION 
3.1 Line profile analysis 

We identified the hydrogen recombination lines, such 
as Ha and H/3, and collisionally excited lines, such as 
[N 11] AA6548,6584, [S 11] AA6717,6731 present in the spectra, 
making use of the known redshift, z=0. 021371, for Haro 15 
(de Vaucouleurs et al. 1991). The strong lines were used 
to analyse the structure of velocity profiles as they allow 
us to verify the existence of more than one component as 
described in Firpo et al. (2010). The adopted laboratory 
wavelengths were taken from the work of Garcia-Rojas et al. 
(2005). 

From the echelle calibrated spectrum, we cut the wave- 
length range where a given emission line is and we trans- 
formed from wavelength to velocity plane using the Doppler 
correction. By measuring the central velocity (wavelength) 
and width of several emission lines we determine the radial 
velocities and velocity dispersions of the ionised gas in the 
differents star- forming regions of Haro 15. The radial veloc- 
ity and the intrinsic velocity dispersion (aint), corrected for 
the instrumental and thermal contributions of each emis- 
sion line, are also derived. For these observations we con- 
sider (Ji = 5.2 kms~^ as instrumental width. The thermal 
contribution was derived from the Boltzmann's equation 
{at = 2kT/ma), where k is the Boltzmann's constant, T the 
kinetic temperature (T ~ 10^ K) and rua the atomic mass of 
the corresponding element. Although a detailed determina- 
tion of the electron density and temperatures, and chemical 
ionic and total abundances for each region will be presented 
in Hagele et al. (2011, in prep., hereafter Paper II), small 
changes of a few hundred degrees do not noticeably mod- 
ify this correction, as previously discussed in Firpo et al. 
(2010). The fluxes were derived from the amplitude (A) and 
the FWHM of the Gaussian profile obtained in the compo- 
nent fitting (F=1.0645x AxFWHM), and the corresponding 
errors were estimated taking into account the errors in these 
two parameters. 

As already reported by Firpo and collaborators, in the 
present work we have also found that all Haro 15 emission 
knots show evidence of wing broadening, which is always 
found in the Ha line profiles and it is usually also observed 
in the profiles of the bright emission lines. Making use of the 
iterative fitting of multiple Gaussian profiles we evaluated 
the presence of a broad component and more than one nar- 
row component present in the emission line profile. In this 
case, we fitted a broad component, explaining the integral 
profile wings for all regions. 

In the following subsections we will discuss our findings 
for each knot resulting from the profile fits. 
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Ham 15 A 

In our high resolution spectra Haro 15 A shows a complex 
structure which, although evident in radial velocity space, 
could not be spatially resolved. 

We identified and fitted Gaussian profiles to the H/3, 
[O iii]A5007, [N ii]A6548, Ha, [N ii] A6584 and [S ii] AA6717 
lines in this knot with ngaussfit routines. The Gaussian 
fits of these emission line profiles are shown in Figures 2 
and 3. The Gaussian fits in the profiles reveal the presence 
of two distinctly separated kinematical components labelled 
narrow 1 and 2 (nl and n2). Component nl shows a pro- 
file slightly broader than component n2, (aint — 28 and 24 
kms~^, respectively) and both components show a relatively 
large spread in individual radial velocities among the differ- 
ent emission lines present in the spectrum. The reliability of 
these values is confirmed when we improve the profile fitting 
using NGAUSSFIT which in turn yields values for the profile 
width of each component. It is worth noting that, although 
the profile fitting to the [O ill] lines shows similar overall 
results, the n2 components show the broader profile, oppo- 
site to what is found for the other emission lines. This could 
be related to a different kinematic behaviour of the highly 
ionised gas, although this needs to be confirmed for other 
high excitation lines. 

Always considering the presence of two distinct compo- 
nents with different radial velocities, the overall fit continues 
showing the presence of a residual in the emission line wings. 
Following the procedures outlined in Firpo et al. (2010), we 
are able to fit a broad component, with a velocity dispersion 
of about 78 kms~^ from the Ha emission line and slightly 
lower from the rest of the lines. Table 2 shows the parame- 
ters for the three components that fit the global profile. In- 
dividual Gaussian component fluxes are listed as fractional 
emission measures (EM/) relative to the total line flux fol- 
lowing the work by Relano & Beckman (2005). The sum of 
these individual fluxes, which we will hereafter refer to as 
overall Ha flux, uncorrected for reddening, is found to be 
(5.49 ± 0.05) X 10~^^ ergs~^cm~^. In Figure 2 we show the 
NGAUSSFIT fitting done with three different Gaussian com- 
ponents in the emission lines which have enough signal to 
provide a reliable fit. The validity of the profile multiplicity 
and broadening is checked over the different emission lines, 
becoming more evident for the strongest emission lines. Fig- 
ure 3 shows the excellent agreement among individual fits 
for the most intense emission lines in Knot A. 

The top panel in Figure 4 shows the spatial profile of the 
Ha emission line in knot A. Two distinct zones are clearly 
distinguished, separated by 1.72 arcsec 0.72 kpc. We have 
performed additional extractions of these two zones (labelled 
A I and A II) individually. Inspection of the individual plots 
displayed in the bottom panels show that the kinematical 
structure is dominated by region A I, somewhat expected as 
it is brighter than region AH, which exhibits a more simple 
emission profile in turn. High spatial resolution integrated 
field spectroscopy is needed to disentangle this complex be- 
haviour. 

Haro 15 B 

The spectrum of region B is dominated by nebular emission 
lines, where we detected and measured [O ill] A4363 and 



He II A4686. We identified and fitted Gaussian profiles to 
the H/3, [O III] AA4959,5007, [N ii] A6548, Ha, [N ii] A6584 
and [S ii] AA6717,6731 lines in this region using ngaussfit. 
It is interesting to note that, although from inspection of 
the Wide Field Planetary Camera 2 image shown in Figure 

1, a complex structure of region B is apparent, this is not 
refiected in the kinematical behaviour of the region, since it 
is possible to fit a single Gaussian component with a veloc- 
ity dispersion of 20 kms~^. Nevertheless, the profile fittings 
show the presence of a residual in the emission line wings 
that we were able to explain by fitting a broad component 
with a velocity dispersion of about 43 kms~^ Results of the 
fitting procedure can be seen in Figure 5; it is worth noting 
that the logarithmic scale chosen for the y-axis magnifies the 
1% residual in the profile wings for Ha and [O ill] A5007. 
For the listed lines the derived radial velocities, velocity dis- 
persions and their corresponding errors are listed in Table 

2. The overall Ha fiux, uncorrected for reddening, is found 
to be (9.88 ± 0.3) x 10"^^ ergs-^cm-^ 

As concluded by Lopez-Sanchez & Esteban (2009), knot 
B has a much higher ionization degree than the rest of the 
regions in Haro 15. These authors suggest that this may be a 
consequence of the extreme youth of this knot. This fact and 
the low metallicity found in knot B (lower than for the other 
regions) indicate that it has a different nature, probably be- 
ing the remnant of a dwarf galaxy which is experiencing a 
minor interaction with Haro 15 (Lopez-Sanchez 2010). 

Haro 15 C 

Although the observed spectrum of region C is faint, we 
were able to identify and fit Gaussian profiles to the H/3, 
[O III] A5007, Ha, [N ii] A6584 and [S ii] AA6717,6731 lines. 
Knot C shows one narrow component together with the un- 
derlying broad component. However, due to the poor signal 
of the spectrum, all emission lines were fitted using the Ha 
NGAUSSFIT solution as a template, fixing the profile centers 
and widths, and only allowing the task to fit the profile am- 
plitudes. In Figure 6 we only show Ha since the rest of the 
intense lines are very noisy due to the low signal- to- noise ra- 
tio of the spectrum. Table 3 shows the derived kinematical 
parameters. The overall Ha flux, uncorrected for reddening 
is (2.67 ± 0.07) X 10"^^ ergs^^cm^^ 

Haro 15 E 

We detect the presence of a broad proflle, suggesting the 
presence of a double peak together with a flux excess in 
the strongest emission lines. We identifled and fltted Gaus- 
sian profiles to the H/3, [O ill] A5007, Ha, [N ii] A6584 and 
[S ii] AA6717,6731 lines. The results of the profile fitting pro- 
cedure yields two components of similar width. Although 
both are supersonic profiles it is not possible to distinguish 
a broad and narrow component and therefore refer to the 
identified components as narrow 1 and 2 (see Figure 6). For 
the listed lines the derived radial velocities, the velocity dis- 
persions and the corresponding errors are listed in Table 3. 
The overall Ha fiux, uncorrected for reddening, is found to 
be (4.99 ± 0.14) x 10"^^ ergs"^cm"^ 

Owing to the poor signal-to-noise spectrum we fit- 
ted some of the emission lines such as [N ii] A6584 and 
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Figure 2. Strongest emission lines from Haro 15 A spectrum. Each panel includes flux calibrated spectrum, where the individual x-axis 
have been normalised to the observed radial velocity for comparison purposes. To enhance details at low luminosity levels, the y axes 
are shown in logarithmic scale. From top to bottom and left to right: H/3, [O ill] A5007, Ha, [N ll] A6584 and [S ll] A6717. 
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[S ii] AA6717,6731, using the Ha ngaussfit solution, and 
iterating separately for individual set of parameters, follow- 
ing the procedure described in Firpo et al. (2010) . 



Haro 15 F 

The observed spectrum of knot F is the weakest one, and 
hence we only identified and fitted the Gaussian profiles of 
the H/3, [O III] A5007, and Ha lines. H/3 and [O iii] A5007 
are in fact very noisy and narrow, rendering impossible the 
task of fitting multiple Gaussian components (narrow plus 
broad). The spectral line broadening in their integrated 
spectra is only observed in the Ha line where we could 
fit a narrow Gaussian component together with the broad 
one, with velocity dispersions of about 8.3 kms~^ and 22 
kms~^, respectively (see Figure 6). The estimated velocity 
dispersion for the narrow component in the Ha line is sub- 



sonic (8 kms~^), typical of classic Hii regions, although the 
presence of a broad component is almost exclusive of Giant 
Hii Regions. 

For Ha, the derived radial velocity, the velocity disper- 
sions and the corresponding errors are listed in Table 3. The 
overall Ha flux, uncorrected for reddening, is found to be 
(1.06 di 0.16) X 10-^^ ergs-^cm-^ 

3.2 Radial Velocities 

The radial velocities for each component of the different 
knots of Haro 15 are given in the corresponding tables where 
the results for the profile fits are shown. In Table 4 we list the 
average radial velocity for each component, together with 
the corresponding errors. The single component is the result 
of a single Gaussian fit to the emission line profiles. In order 
to compare our results with those of Lopez-Sanchez & Es- 
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Figure 3. Same set of lines plotted in Figure 2 now overlapped in the same diagram. To enhance details at low luminosity levels, the y 
axes are shown in logarithmic scale. The evident similarity among different lines belonging to different elements confirms the presence 
of a complex kinematical structure. The results of fitting of individual components for the [S ll] A6717 forbidden line is included as a 
reference. 



1 

Haro 15 A 


1 ' 1 ' 

- HP 

[O m] %5QQ1 
"""^^ - Ha 

\ [N n] A6584 

- [Sn]A^717 








1.1,1. 



"< 10 



Velocity (km s ) 



Figure 4. Top panel: Intensity distribution along the slit at the peak of the Ha emission. Bottom panels: Result of profile extraction at 
the locations of regions A I and A II, together with the result of individual Gaussian profile fitting. To enhance details at low luminosity 
levels, the y axes are shown in logarithmic scale. Two narrow components and a broad one are fitted for A I and only a single narrow 
and broad component are identified in A II. 
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teban (2009) we have plotted our radial velocities in Figure 
7 relative to the velocity of the galactic center 6415 kms~^ 
(Lopez-Sanchez 2006), after correcting to the Galactic Stan- 
dard of Rest (GSR) using the rvcorrect task by IRAF. 
Figure 7 shows our results together with data presented by 



Lopez-Sanchez & Esteban (2009) for their long-slits with 
position angles (PA) of 41° and 117°. Regarding the single 
component velocity, it can be clearly seen that the aver- 
age radial velocities derived from the emission lines follow 
nicely those derived by Lopez-Sanchez & Esteban (2009). 
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Figure 5. Ngaussfit fits with two Gaussian components in the Haro 15 B strongest emission hne profiles, a narrow component and 
a broad one. In order: H/3, [O ill] A5007, Ha, [N ll] A6584 and [S ll] A6717. Note that the logarithmic y-scale magnifies the actual low 
(1%) error in the overall fits 
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However, it is worth mentioning that individual velocities of 
narrow components (nl and n2) for knot A differ by about 
60 kms~^. This cannot be easily explained by scatter within 
galactic rotation, suggesting that their relative velocities are 
due to mutual orbital motion around their gravity center. 



3.3 Relation between Ha Luminosities and 
Velocity Dispersion 

Figure 8 shows the distribution of Haro 15 knots in the 
L(Ha) vs. a plane. Each knot is labelled from A to F as 
used throughout this paper. Narrow components are identi- 
fied with subscript n while subscript s refers to line widths 
measured by fitting a single Gaussian component to the line. 
Also plotted are data on NGC6070 and NGC 7479 from 
Firpo et al. (2010) and Giant Hii Regions from Bosch et al. 
(2002) together with the regression line fitted by those au- 



thors to their "young" Giant Hii Regions. Ha luminosities 
were derived from the fluxes measured directly from the 
component fitting to our echelle spectra, and using distances 
as published by de Vaucouleurs et al. (1991) for Haro 15. The 
correction for reddening could change flux values by a fac- 
tor up to 1.72 (0.24 in logarithm) (Lopez-Sanchez & Esteban 
2009) but, as we cannot determine it for every regions, we 
prefer to list the uncorrected ones. 

Similar to what was found by Firpo et al. (2010), the lo- 
cation of the studied regions is far from random. Haro 15 Hii 
regions show a correlation between their luminosities an ve- 
locity dispersions. If we focus on the fit of single components 
they seem to follow the regression found for virialised sys- 
tems. The outstanding exception to this trend is region A, 
but this can be attributed to the fact that the profile can be 
split in two narrow components, resembling the behaviour 
found for regions NGC 7479-1 and NGC 7479-H (Firpo et al. 
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Table 2. Results of Gaussian profiles fitting to the observed emission lines in Haro 15 A and B. Each emission line is identified by its 
ion laboratory wavelength and ion name in columns 1 and 2. According to the different fits performed on each line, column 3 identifies 
each "narrow component (1 and 2, where applicable), and a broad component. Radial velocities (Vr) and intrinsic velocity dispersions 
(o'int) together with their respective errors are expressed in kms"-*^. The intrinsic velocity dispersions are corrected for the instrumental 
and thermal widths. Fractional emission measures (EMf) in %. 



Ao 


ion 


comp. 






Haro 15 
error aint 


A 

error 


EM/ 




Haro 15 
error aint 


B 

error 


EM/ 


4861 


H/3 


narrow 


1 


6375.5 


0.8 


28.8 


1.1 


55.7 


6381.2 


0.8 


17.3 


1.1 


23.6 






narrow 


2 


6440.1 


0.9 


20.7 


1.1 


25.6 
















broad 




6390.6 


5.2 


72.5 


2.4 


18.7 


6393.4 


0.8 


43.3 


1.1 


76.4 


4959 


[O III] 


narrow 


1 


6373.4 


1.1 


26.5 


1.3 


39.0 


6383.2 


0.8 


16.5 


1.1 


15.5 






narrow 


2 


6426.2 


1.3 


31.9 


1.5 


38.0 
















broad 




6384.4 


4.7 


64.0 


6.0 


23.0 


6394.4 


0.7 


40.1 


1.1 


84.5 


5007 


[O III] 


narrow 


1 


6372.7 


0.7 


24.9 


1.1 


39.2 


6384.1 


0.8 


15.1 


1.1 


13.0 






narrow 


2 


6429.1 


0.8 


28.8 


1.1 


35.4 
















broad 




6396.2 


1.6 


73.0 


2.4 


25.4 


6395.2 


0.7 


40.2 


1.0 


87.0 


6548 


[N II] 


narrow 


1 


6378.7 


1.2 


28.6 


1.1 


53.4 


6381.9 


1.3 


20.6 


1.2 


50.2 






narrow 


2 


6438.3 


1.1 


21.5 


1.3 


27.9 
















broad 




6405.5 


12.1 


72.0 


14.1 


18.7 


6395.6 


4.2 


49.0 


4.7 


49.8 


6563 


Ha 


narrow 


1 


6376.1 


0.6 


27.3 


0.8 


47.6 


6385.0 


0.6 


20.0 


0.8 


28.4 






narrow 


2 


6436.0 


0.6 


24.5 


0.8 


33.0 
















broad 




6395.8 


2.6 


77.8 


1.7 


19.4 


6393.5 


0.6 


43.7 


0.8 


71.6 


6584 


[N II] 


narrow 


1 


6379.9 


0.6 


27.5 


0.8 


46.6 


6381.5 


0.7 


24.3 


0.9 


59.8 






narrow 


2 


6441.9 


0.7 


24.5 


0.9 


28.1 
















broad 




6399.2 


2.1 


74.7 


2.5 


25.3 


6392.9 


2.3 


48.4 


2.3 


40.2 


6717 


[S II] 


narrow 


1 


6379.1 


1.2 


31.7 


0.9 


38.1 


6382.0 


0.8 


25.1 


0.8 


67.1 






narrow 


2 


6441.8 


1.2 


19.1 


1.3 


16.5 
















broad 




6397.9 


2.2 


65.5 


2.4 


45.5 


6387.8 


3.2 


52.4 


3.3 


32.9 


6731 


[S II] 


narrow 
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6378.3 


1.5 


32.3 


1.7 


37.2 


6384.9 


0.8 


23.4 


1.0 


66.6 






narrow 
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6442.3 


1.1 


16.7 
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20.2 
















broad 




6398.3 


3.4 


65.3 


4.3 


42.6 


6380.6 


4.8 


59.0 


5.3 


33.4 



2010). The narrow components of the composite fits show, as 
expected, relatively smaller luminosities and velocity widths, 
but they still lie around the same linear regression. In this 
case, it is region F that fails to follow the trend but, as men- 
tioned earlier, its low velocity dispersion makes it a very 
unlikely candidate for a Giant Hii Region. It can be seen 
from the figure that while narrow components tend to clus- 
ter around the L(Ha) oc relation, single component fits 
would provide a flatter power- law exponent. This is proba- 
bly due to the broad component contributing a substantial 
part of the total Ha flux. Although not included in the fig- 
ure, the broad components are located, as expected, in a 
parallel sequence shifted towards lower luminosities and/or 
higher velocity dispersions. 



4 SUMMARY 

From high resolution spectroscopy of Haro 15 we have per- 
formed a thorough analysis of the emission line profiles of 
several knots considering multiple component fits to their 
profiles. Our results can be summarised as: 

• Giant Hii Regions of Haro 15 show a complex structure 



within the profile of all their emission lines, detected both 
in recombination and forbidden lines. 

• The brightest emission lines can be split in at least two 
strong narrow components plus an underlying broad com- 
ponent. 

• Although regions tend to follow the galaxy kinematics, 
the narrow components of knot A have relative velocities 
that are too large to be explained by galactic rotation. This 
behaviour can be explained if we are observing the orbital 
motion around their gravity center. 

• Almost all knots follow the relation found between lumi- 
nosities and velocity dispersions for virialised systems, either 
when considering single profile fitting or the strong narrow 
components in more complex fits. Among these, the single 
one shows a relatively flatter slope. 

The presence of more that one component in the Gaus- 
sian fits to the emission line profiles, such as those anal- 
ysed in this paper, has been discussed in several previous 
studies. Hagele et aL, (2007; 2009; 2010) showed fits that 
involved the existence of broad and narrow components for 
the emission lines of the ionised gas in circumnuclear star- 
forming regions. Furthermore, Firpo et al., (2010) were able 
to detect two distinct components within the narrow fea- 
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Table 3. Results of Gaussian profiles fitting to the observed emission lines in Haro 15 C, E, and F. Each emission line is identified by 
its ion laboratory wavelength and ion name in columns 1 and 2. According to the different fits performed on each line, column 3 identifies 
the narrow and broad components. Radial velocities (Vr) and intrinsic velocity dispersions ((Tint) together with their respective errors 
are expressed in kms~^. The intrinsic velocity dispersions are corrected for the instrumental and thermal widths. Fractional emission 
measures (EMj) in %. 
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error 
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Haro 15 
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Ha 


narrow 




6494.2 


0.4 


11.8 


0.9 


44.5 






broad 




6492.7 


0.7 


28.5 


1.0 


55.5 










Haro 15 


E 








4861 




narrow 


1 


6303.0 


0.8 


21.4 


1.2 


48.6 






narrow 


2 


6314.4 


1.1 


28.2 


1.4 


51.4 


5007 


[O III] 


narrow 


1 


6302.9 


1.7 


23.9 


2.2 


32.8 






narrow 


2 


6307.6 


1.2 


27.6 


1.5 


67.2 


6563 


Ra 


narrow 


1 


6303.0 


0.7 


21.4 


0.9 


59.1 






narrow 


2 


6314.4 


1.0 


28.2 


1.2 


40.9 


6584 


[N II] 


narrow 


1 


6306.3 


2.0 


22.4 


2.9 


31.2 






narrow 


2 


6312.4 


2.0 


34.2 


2.0 


68.8 


6717 


[S II] 


narrow 


1 


6309.0 


1.9 


17.6 


1.6 


24.9 






narrow 


2 


6314.3 


1.4 


27.8 


1.5 


75.1 


6731 


[S II] 


narrow 


1 


6306.1 


3.5 


14.9 


2.3 


29.6 






narrow 


2 


6314.7 


3.4 


30.0 


2.9 


70.4 










Haro 15 


F 








6563 


Ha 


narrow 




6497.7 


1.1 


8.3 


1.0 


75.9 






broad 




6479.1 


3.4 


22.0 


5.7 


24.1 



Table 4. Average LSR radial velocities for Haro 15 emission line knots. Column 1 indicates the fitted feature and the remaining columns 
list the average radial velocity for each knot and its uncertainty, both in kms"-*^. 

Comp. < Vr > error < Vr > error < Vr > error < Vr > error < Vr > error 

Haro 15 A Haro 15 B Haro 15 C Haro 15 E Haro 15 F 

single 6420.3 0.7 6410.2 0.6 6514.4 0.6 6327.5 0.6 6516.3 0.7 

narrow 1 6397.3 0.9 6403.5 0.8 6516.5 1.3 6325.4 1.4 6518.3 1.2 

narrow 2 6457.1 0.9 ... ... ... ... 6333.2 1.3 

broad 6418.4 4.4 6414.1 1.8 6511.5 4.8 ... ... 6499.8 3.5 



Table 5. Intrinsic velocity dispersions (dint) and the Ho; luminosity (L(¥La)s) for a fitted single Gaussian component, together with 
the their associated errors, expressed in kms"-*^ and ergs"-*^, respectively. We add the Ha luminosity (L(}ia) oyerall) obtained from the 
overall Ha flux, uncorrected for reddening, and its error, in ergs"-*^. 



knot aint L(Ha)s L{'Ha)oyerall 

Haro 15 A 47.3±0.9 4.9±0.7 x 10^^ 4.9±0.7 x 10^^ 
Haro 15 B 34.4±0.8 8.6±1.2 x 10^° 8.9±1.3 x 10^° 
Haro 15 C 19.3±0.8 2.3±0.3 x 10^9 2.4±0.3 x 10^9 
Haro 15 E 24.4±0.9 4.5±0.6 x 10^9 4.5±0.6 x 10^9 
Haro 15 F 11.6±0.8 9.1±1.3 x 10^^ 9.5±1.3 x 10^8 
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Figure 7. Position vs. Velocity diagrams for the observed knots in Haro 15. Ho; radial velocity results, relative to knot A, are over 
imposed on the data by Lopez-Sanchez &: Esteban (2009) for a slit positions PA 41° (left panel) and PA 117° (right panel). The Ho; and 
[O III] 5007 data of Lopez-Sanchez & Esteban (2009) in Haro 15 are identified by blue circles and red triangles, respectively. 
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ture of the emission lines in the brightest Giant Hii Regions 
NGC 7479-1 and IL Although we still lack the number of ob- 
jets needed to support a generalization of this trend, these 
data show that the presence of multiple kinematical com- 
ponents among extragalactic star-forming regions cannot be 
overlooked. They might have a strong impact on subsequent 
analysis that rely on basic parameters, such as the velocity 
dispersion and chemical abundance of the ionised gas, the 
inferences about the nature and strength of the source of 
ionization, or the classification of the activity in the central 
regions of galaxies. 



It is necessary to map these regions with high spectral 
and spatial resolution and much better S/N ratio to dis- 
entangle the origin of these different components. 2D spec- 
troscopy performed with Integral Field Unit (IFU) is the 
ideal tool to tackle this issue. 
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Figure 8. log(L) — log(cr) relation for our HII regions. Luminosities and velocity dispersions are derived from our spectrophotometric 
data. The plot includes results from narrow components (labelled as n) and single profile fits (labelled as s). We add the data from Firpo 
et al. (2010) (n, nA and nB for narrow, and g for single components) identified by color solid error bars and with numbers (from 1 to 
6): NGC 7479-1 (1) in red, NGC 7479-11 (2) in green, NGC 7479-III (3) in yellow, NGC 6070-1 (4) in maroon, NGC 6070-11 (5) in violet 
and NGC 6070- IV (6) in magenta. A few Giant Hll Regions from Bosch et al. (2002) (blue dashed error bars) together with their linear 
fit to their "young" Giant Hll Regions are plotted as a reference value. The luminosities are not corrected for extinction. 
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